Resistivity model in both paramagnetic and ferromagnetic phases 
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The resistivity model as a function of temperature and ionization energy (doping) is derived 
with further confinements from spin-disorder scattering in ferromagnetic phase. Magnetization and 
polaronic effects capture the mechanism of both spin independent and spin-assisted charge transport 
of ferromagnets. The computed Tcrossover below Tc and carrier density in Gai_2,Mna;As system are 
8-12 K and 10^^ cm"'', remarkably identical with the experimental values of 10-12 K and 10^*-10^° 
cm"'' respectively. 
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1. INTRODUCTION 

Diluted magnetic semiconductors (DMS) have the 
tremendous potential for the development of spintron- 
ics and subsequently will lay the foundation to realize 
quantum computing. This applicability arises due to fer- 
romagnetic nature of DMS. In other words, both the 
charge and spin of the electrons can be exploited with 
limited Mn doping in GaAs semiconductor. In order 
to achieve this, one needs to understand the transport 
mechanism such as the variation of resistivity with tem- 
perature and doping in both above and below Tc consis- 
tently. Interestingly, Van Esch et al. have proposed 
multiple exchange interactions, which are ferromagnetic 
(FM) hole-hole and antiferromagnetic (AFM) Mn-hole 
interactions for DMS. These two effects, after neglecting 
the direct exchange between Mn-Mn (due to very diluted 
nature of DMS) are seem to be sufficient enough to de- 
scribe the temperature dependent magnetization curves 
(M(T)) accurately. However, even after inclusion of FM 
and AFM effects including the spin disorder scattering, 
the transport property in the FM phase is still not well 
understood. Unfortunately, this is also true for the case 
of metallic property below Tc in the well known and ex- 
tensively studied FM manganites as pointed out by Ma- 
hendiran et al. 2\. The resistivity {p(T)) above Tc for 
manganites is found to be in an activated form described 
by the equation [2], 



p{T > Tc) ^ po exp 



Eg 

knT 



(1) 



Ea is the activation energy, po and ks denote the resid- 
ual resistivity sXT ^ Ea and Boltzmann constant respec- 
tively. In the FM phase, the influence of M{T)/Mi,2 is 
more pronounced than the electron-phonon {e-ph) contri- 
bution where the latter requires an overwhelmingly large 
coupling constant Note that M4.2 is the magnitude 
of magnetization at 4.2 K. Therefore, Mahendiran et al. 
have suggested that conventional mechanism namely, e- 



ph scattering has to be put aside so as to explain the p{T) 
for manganites below Tc. On the contrary, p{T) with e- 
ph involvement for DMS in the paramagnetic phase is 
given by 

^(j., r )- C,+C,[e^p{en/T)-l]-' 

^ ^' kBT\u[l + e^^{{Em-Ef)/kBT)y ^ ' 

The term, C2/ [exp(8j3/T) — l] takes care of the e- 
ph contribution. O^, Ef, Ci and C2 represent 

the Debye temperature, Fermi level, mobility edge and 
numerical constants respectively. The p{T) in the FM 
phase based on the spin disorder scattering as derived by 
Tinbergen-Dekker is given by 0| 



PsoiT < Tc) = 



(m;,)5/2^(2j;^)l/2 
7r(n,p)e^7i^ 



S{S + \)-S- 



21 Mtd{T) 



Ma 



' ) [2TSiS+l)M,. 



(3) 



TV is the concentration of nearest neighbor ions (Mn's 
concentration) while {n,p) is the concentration of charge 
carriers (electrons or holes respectively), m* ^ denotes ef- 
fective mass of electrons or holes, h = h/2'n, h = Planck 
constant, e is the charge of an electron, Ep and Jex are 
the Fermi and FM exchange interaction energies respec- 
tively while S is the spin quantum number. Equation (PJ 
becomes equivalent to Kasuya 4] if one replaces the term, 
teinh[3TcMTD{T)/2TS{S +1)M4,,2] with 1. Again, an 
accurate equation for the p{T) below Tc is still lack- 
ing since spin disorder scattering alone is insufficient as 
shown bv Tinbergen and Dekker Q as well as reviewed by 
Ohno p|. Hence, it is desirable to derive a formula that 
could describe the transport mechanism of ferromagnets 
for the whole temperature range i.e., for both paramag- 
netic and FM phases and even at very low T. With this in 
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mind, the ionization energybased Fermi-Dirac statistics 
(iFDS) as derived in Ref. jl S E EH and spin dis- 
order scattering based resistivity models derived by both 
Tinbergen-Dekker and Kasuya will be employed in order 
to derive p as a function of T, ionization energy (Ej) and 
Mp(r, M4.2). The consequences of p{T, Ej, Mp{T, M4.2)) 
that arises from the variation of T, Ej and Mp{T, M4.2) 
are discussed in detail based on the experimental data 
reported by Van Esch et al. Q and Mahendiran et al. ■ 

2. RESISTIVITY MODEL 



measures the combination of electrons and its strain 
field due to neighboring ions, which is nothing but po- 
larons 0, ■ The absolute value of Ej can be obtained 
from Pj -E'/ = e^/STreeors. e and eo are the dielectric 
constant and permittivity of free space respectively, vb 
is the Bohr radius. Furthermore, the variation oi Ej with 
magnetic field, H will give rise to an inverse variation on 
vb that also takes care of the polaronic effect 0, 
Substituting I/ts — AT^ (due to electron-electron inter- 
action), Eqs. and (0) or © into Eq. Q), then one can 
arrive at 



The total current in semiconducting ferromagnets with 
contributions from both paramagnetic and FM phases is 



z\hK sh\ For convenience, the 



spin-up, t denotes the direction of the magnetic field or a 
particular direction below Tc, while the spin-down, J, rep- 
resents any other directions. Note that the total energy 
(Kinetic -I- Magnetic), -Sj^+jv/ ^k+m energy 
level splitting below Tc- As such, the total current can 
be simplified as J = -f Jj^ = Jg + Jse if the considered 
system is a n-type while J = Jh + Jsh if it is a p-type. 
Je and Jh are the spin independent charge current (elec- 
trons and holes respectively) in the paramagnetic phase 
whereas Jse and Jsh are the spin-assisted charge current 
in the FM phase. Thus the total resistivity (n or p-type) 
can be written as 



Pe,h + Pse.sh 



K.h 


-1 

+ 








{n,p)e^TsD_ 



(4) 



TsD represents the spin disorder scattering rate. The 
carrier density for the electrons and holes (ri,p) based on 
iFDS are given by 0,1,11 [H El 



3/2 



(m*)^/^ exp 



Ep — Ej 



3/2 



{fUhf/'^ exp 



-Ep — El 
knT 



(5) 



(6) 



The derivation of iFDS, f{Ei) — exp [ — /i — 
HEinitiai state ±Ei)] by employing the restrictive Condi- 
tions, dUi = and Y.T {^initial state ± Ei)idni = 

is well documented (including its applications) in the 
Refs. imilHIHinillllll. Einitiai state dcuotcs the 
energy at certain initial state and the Lagrange multipli- 
ers, He + XEi = - In [{n/V){2TrXh^/mef/^], Hh^^^Ei = 
In [{p/V){2T:Xn^ /ruhf'^] and A = l/fcsT. V is the vol- 
ume in k space and the gap-parameter, Ei that rep- 
resents electron-ion attraction is also a parameter that 
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e^sei'E) — 

which, 
= 2m* 



ABe^p [{Ei + EF)/kBT] 
AT^/^[Mp{T, M4.2)]-i + BT-1/2 ■ 



(7) 



A 



= [AeM/2e^ml^y/^][2nh'/kBr^^ 



^NinEpf'^JlJe'nkT and 



' SD 



B 

[N{2EFf'''{mlhf/''/TTnyi^Mp{T,Mi.2). Ae,h is the 
T independent electron-electron scattering rate constant. 
The empirical function of the normalized magnetization 
is given by 



MpiT,AU,2) = 1 



MpjT) 



M. 



4.2 



(8) 



Equation (jS)) is an empirical function that directly 
quantifies the influence of spin alignments in the FM 
phase on the transport properties of charge and spin 
carriers 0| in accordance with Eq. (jjj). In other 

Al (T) 

words, the only way to obtain j^'^^ ^ is through 
Eq. IHl. In fact, Eq. |(Hl is used to calculate 
MTD{T)/Mi,2 and Mk{T) /A'U,2 by writing S{S + 

' Mtd{T) \2 _ _5.|-MTD(T)U„„k r 3TcMtd(T} 



Ma 



Ma 



tanh [ 



2TS(S+l)Mi,2- 



Mp{T, M4.2) and S{S + 1) - 5^ 



Lpy-L,iv^i.-2) o.,,^ uyu ^ O ( jv/4.2 ) ' ^ \ M4.2 > " 

Mp{T, M4.2) respectively. Consequently, one can actually 
compare and analyze the Ma{T)/M4,2 (a = TD, K, p) 
calculated from Tinbergen-Dekker (TD), Kasuya (K) and 
Eq. with the experimentally measured Mexp{T)/Mi,2- 
However, one has to switch to Eq. © given below for the 
hole-doped strongly correlated non-ferromagnetic semi- 
conductors, which is again based on iFDS U0i 



Ph 



Ah{ml 



2e2 



3/2 



exp 



El -\- Ep 



knT 



(9) 



3. DISCUSSION 
3.1. Diluted magnetic semiconductors 

It is the purpose of this paper to explain the trans- 
port mechanism below Tc without violating the physi- 
cal properties known in p(T > Tc). The discussion on 



3 



the mechanism of transport properties of ferromagnets is 
seen through the eyes of resistivity measurements both 
in the presence of and in the absence of H. The resistiv- 
ity measurements [J and its fittings based on Eqs. Q) 
and are shown in Fig. a) and b) respectively for 
Gai-ajMuj^As. Literally, one needs two fitting param- 
eters {A and Ej) for p{T > Tc) and another two {B 
and Mp(r,M4.2)) for p(T < Tc). All the fitting param- 
eters are given in Table Note that 5=1 and 5/2 
are employed for the fittings of Mk(T)/M4,2 while Tc 
and Tcr were determined from the experimental resis- 
tivity curves, not from the magnetization measurements 
or any other techniques. The deviation of Mj^ (T)/M4.2 
from the Mf.xp/M4 2 increases with S from 1 — > 5/2. The 
p{T) is found to increase with x from 0.060 to 0.070 due 
to the mechanism proposed by Van Esch et al. ^^'^ 
Ando et al. 0|. They proposed that neutral Mn^+ ac- 
ceptors that contribute to magnetic properties could be 
compensated by As, where for a higher concentration of 
Mn, instead of replacing Ga it will form a six-fold coordi- 
nated centers with As (Mn^-^^) [lElEl- These centers 
will eventually reduce the magnitude of ferromagnetism 
(FM) in DMS due to the loss of spin-spin interaction be- 
tween Mn(3d^) and h. Parallel to this, Mn^'^" formation 
is substantial in such a way that Mn'^"'' ions do not sub- 
stitute Ga^+ ions. Therefore, p{T) will be influenced by 
]y[j^6As clusters, defects and Ga-Mn-As phase simultane- 
ously significantly in this range of x. This is also indeed 
in fact in accordance with iFDS based resistivity mod- 
els since if one assumes Mn^+ {Ej — 1113 kJmol^^) or 
Mn3+ {Ei = 1825 kJmol" i) substitutes GaP+ {Ej = 1840 
kJmol"^), then p(T) should further decrease Q with x, 
which is not the case here. Thus, iFDS also suggests 
that Mn^+ or Mn'^+ do not substitute Ga^+. Interest- 
ingly, the Tc-s observed in Gao.94oMno.o6oAs {T^ = 10 
K, annealed: 370°C) and Gao.gaoMno.oToAs (r„. = 12 K, 
as grown) are identical with the calculated values, where 
El + Ep = 8 K and 12 K respectively. Ei + Ep is actu- 
ally equivalent to Tcr because of its exponential contribu- 
tion as shown in Eq. The calculated carrier density 
using Ej + Ep (8, 12 K), mjjj — rest mass and Eq. © 
is 2.4 X 10^^ cm^"^. Below Tc, spin alignments enhance 
the contribution from J^e and reduces the exponential 
increase of p{T). This reduction in p{T) is as a result of 
dominating Jge and small magnitude of Ej + Ep (8-12 
K), consequently its effect only comes around at low T as 
clearly shown in Fig. ^a). The Gap. 930 Muq. 070 As sam- 
ples after annealing at 370 °C and 390 °C do not indicate 
any FM ^ (Fig.[T]b)). Thus the fittings are carried out 
with Eq. (jSJ that only require two parameters namely, A 
and Ej + Ep since Jge — (there is no observable Tc) 
and/or dM^iT) / Mi,2dT = {Mp{T,Mi,2) = constant). 
The exponential increase of p{T) is due to Ei + Ep as 
given in Eq. Q with zilch Jse contribution. 

Figure ^c) and d) indicate the calculated normalized 
magnetization, Ma{T)/M4,2 obtained from Eq. {Tjl. Note 



that Mp^TD,K{T)/Mi,2 is a fitting parameter that has 
been varied accordingly to fit p{T < Tc). As a matter of 
fact, Mp{T,Mi,2) is used to calculate Mp^td,k{T)/Ma,2 
with S = \. Mp,TD.K{T)/Mi,2_\s also compared with 
the experimentally determined [ij Mexp{T)/M4,2 as de- 
picted in Fig. n d). One can easily notice the re- 
lation, MTD{T)/Mi,2 > MK{T)/Mi,2 > Mp{T)/M4,2 
> Mexp{T)/M4,2 from Fig. [Tl c) and d). As such, 
Mp(T)/M4,2 determined from Eq. Q is the best fit for 
the experimentally measured Mexp{T) / M4,2. Obviously, 
the higher number of aligned spins as calculated from 
Eq. using Mp{T)/Mi,2 compared to Mexp{T)/Mi,2 
is due to the ability of both Je and Jse to follow the 
easiest path. Simply put, resistivity measures only the 
path with relatively lowest Ei and with easily aligned 
spins that complies with the principle of least action. In 
contrast, magnetization measurement quantifies the av- 
erage of all the spins' alignments. On the other hand, the 
discrepancies of Mtd{T)/M4,2 with Mexp{T)/M4,2 and 
Mk{T)/M4,2 with Mexp{T)/M4,2 are due to long range 
FM hole-hole and AFM Mn-hole interactions [jj] apart 
from the ability of both Je and Jse to follow the eas- 
iest path. The violation between Mtd.k{T)/M4,2 and 
Mexp{T)/M4,2 suggests that the spin disorder scattering 
alone is inadequate, in which, the principle of least action 
have had played an enormous role. 

3.2. Manganites 

Now switching to manganites, Mahendiran et al. 
discussed p{T < Tc) with respect to Eq. and obtained 
the activation energy, Ea = 0.16 eV for x = 0.1 and 0.2 of 
Lai_j;Gaa;Mn03 samples at T. Using Eq. O however, 
Ej + Ep for the former and latter samples are calculated 
to be 0.12 and 0.11 eV respectively. The calculated car- 
rier density using Ej + Ep (0.12, 0.11 eV), = rest 
mass and Eq. ^ is approximately 10^^ cm"-^. In the 
presence ofH. — 6 T, Ej + Ep is computed as 0.0776 eV 
for X = 0.2 that subsequently leads to p = 10^^ cm""^. It 
is proposed that the activated behavior for p{T > Tc) is 
due to Ej , Coulomb interaction between ion and electron 
or rather due to the polaronic effect IEE3- The fittings 
are shown in Fig.|21a) and b) while its fitting parameters 
are listed in Table HI Theoretically |E[I3> Ca2+ (Ej = 
868 kJmol-i) < La3+ {Ej = 1152 kJmol"!), therefore 
p{T) is expected to decrease with Ca^"*" doping signifi- 
cantly. Gontradicting to that, only a small difference of 
Ei + Ep between x = 0.1 (0.12 eV) and 0.2 (0.11 eV) 
is observed due to Mn'*+ 's compensation effect where the 
quantity of Mn^+ increased 6% from x — 0.1 (19%) to 
0.2 (25%) 0. To clearly see this, the difference of Ej be- 
tween Ca^"*" and La'^"'" is calculated, which is 1152 — 868 
= 284 kJmol"^ and subsequently it is compared with the 
6% increment of Mn3+^4+ ^ 494Q kjniol-i), which 
is 0.81(1825) -I- 0.19(4940) - 0.75(1825) - 0.25(4940) = 
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187 kJmol ^. Consequently, the actual difference is only 
284 - 187 = 97 kJmol-^ instead of 284 kJmor^ This 
simple calculation exposes that Ca^+'s contribution has 
been compensated with 6% additional Mn^+ . All the val- 
ues of Ej discussed above were averaged in accordance 
with Ei[X^+] = J2i=i ^iid should not be taken liter- 
ally since those £^/s are not absolute values. The abso- 
lute values need to be obtained from the rs dependent 
Ei equation stated earlier. Prior to averaging, the 1**, 
2"^^, 3'''' and 4*'' ionization energies for all the elements 
mentioned above were taken from Ref. [l^ . 

As is well known, at 6 T, Lao.8Cao.2Mn03 indicate 
a much lower resistivity (Fig. 13 b)). The result that 
larger H giving rise to conductivity at T > Tc is due 
to relatively large amount of aligned spins at higher T 
or H gives rise to Jse at a higher T. Hence, Tc at 6 
T > Tc at T and one can also conclude, rs (at 6 T) 
> rs (at T) due to the inequality, Ej + Ep = 78 
meV (at 6 T) < Ej + Ep = 112 meV (at T) com- 
plying with Eq. Q and iFDS. Figure [21c) and d) depict 
calculated M„(r)/M4.2 with 5 = 1 and M^^p{T)/Mi,2 
for X = 0.2. Calculated Mtd{T)/M4,2 is dropped for 
Lai_a;Ca^Mn03 since Mk{T)/M4,2 seems to be a bet- 
ter approximation than MrnCT) /M4__2 as indicated in 
Fig. He) and d). Millis et al. [l7| have shown theoret- 
ically that double exchange mechanism (DEM) alone is 
inadequate and one needs to incorporate polaronic effect 
into DEM. In fact, Eq. {T)) takes both polaronic effect 
(Ei) and DEM (M„(T)/M4.2) into account and yet there 
is a discrepancy between Mp{T)/Mi,2 and Mexp{T)/M4,2 
though Eq. O reproduces p{T) at all T range accurately. 
Again, this incompatibility is due to the principle of least 
action as stated earlier. Note that the discrepancy be- 
tween Mk{T)/M4,2 and Mn(T)/MA2 indicate the inad- 
equacy pointed out by Millis et al. [13, not the former. 
This inadequacy may not be clear from Fig.|2]d) because 
the deviation is large with respect to magnetization in 
a narrow range of T (careful observation will reveal this 
though). In addition, the manganites' charge transport 
mechanism below Tc is also in accordance with Eq. (|7|) 
because the term, Mp{T, M4.2) handles the exchange in- 
teractions' complexities separately for DMS and man- 
ganites. For example, one can clearly notice the different 
type of discrepancies between DMS and manganites by 
comparing the empirical function, Ma{T) /Mi,2 (a = P, 
exp) between Fig. [Hd) and Fig. Old). Hence, Eq. is 
suitable for both types of ferromagnets, be it diluted or 
concentrated. In conclusion, the resistivity model that 
have incorporated the polaronic and magnetization ef- 
fects has been derived based on iFDS and spin disor- 
der scattering theories. This model is able to explain 
the transport mechanism below the Curie temperature 
without violating the physical properties above Tc- The 
discrepancy of the magnetization curves calculated from 
this resistivity model with the experimental data of DMS 
arises as a result of the ability of both spin and charge 



TABLE I: Calculated values of T independent electron- 
electron scattering rate constant [A), B, which is a func- 
tion of T independent spin disorder scattering rate constant 
and spin exchange energy (Jex) as well as the ionization en- 
ergy (Ei). Note that Tc and Tcr were obtained from the 
resistivity curves. All these parameters are for Mn doped 
Gai_j;Mn2,As (as grown and annealed at 370 °C, 390 °C) and 
Ca doped Lai_2,Caa;Mn03 (measured at and 6 T) systems. 
All Gai_a;Mna;As samples were measured at T. 



FIG. 1: Equation Q has been employed to fit the experi- 
mental p(r) plots for Gai_a;Mna;As as given in a) whereas 
Eq. is used to fit the plots in b). All fittings are indicated 
with solid lines, b) is actually for annealed non-ferromagnetic 
Gao.93oMno.o70 As samples, c) and d) show the T variation of 
calculated Ma{T)/M4,2 (a = K, TD, p) with S = 1 for a; = 
0.060 and 0.070 respectively. Mk{T)/Ma.2 is also calculated 
with S — 5/2 to indicate the type of deviation one should 
expect with increasing S. The experimental Mexp(T) / M4„2 
plots is only available for x — 0.070 (as grown) as shown in 
d). 

currents to follow the easiest path according to the prin- 
ciple of least action, long range ferromagnetic and short 
range antiferromagnetic interactions. Whereas for man- 
ganites, this discrepancy is still noticeable solely because 
of the principle of least action. 
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